Undiagnosed RASopathies in infertile men.

RASopathies are syndromes caused by congenital defects in the Ras/mitogen-activated protein kinase (MAPK) pathway genes, with a population prevalence of 1 in 1,000. Patients are typically identified in childhood based on diverse characteristic features, including cryptorchidism (CR) in >50% of affected men. As CR predisposes to spermatogenic failure (SPGF; total sperm count per ejaculate 0-39 million), we hypothesized that men seeking infertility management include cases with undiagnosed RASopathies. Likely pathogenic or pathogenic (LP/P) variants in 22 RASopathy-linked genes were screened in 521 idiopathic SPGF patients (including 155 CR cases) and 323 normozoospermic controls using exome sequencing. All 844 men were recruited to the ESTonian ANDrology (ESTAND) cohort and underwent identical andrological phenotyping. RASopathy-specific variant interpretation guidelines were used for pathogenicity assessment. LP/P variants were identified in PTPN11 (two), SOS1 (three), SOS2 (one), LZTR1 (one), SPRED1 (one), NF1 (one), and MAP2K1 (one). The findings affected six of 155 cases with CR and SPGF, three of 366 men with SPGF only, and one (of 323) normozoospermic subfertile man. The subgroup "CR and SPGF" had over 13-fold enrichment of findings compared to controls (3.9% vs. 0.3%; Fisher's exact test, p = 5.5 × 10-3). All ESTAND subjects with LP/P variants in the Ras/MAPK pathway genes presented congenital genitourinary anomalies, skeletal and joint conditions, and other RASopathy-linked health concerns. Rare forms of malignancies (schwannomatosis and pancreatic and testicular cancer) were reported on four occasions. The Genetics of Male Infertility Initiative (GEMINI) cohort (1,416 SPGF cases and 317 fertile men) was used to validate the outcome. LP/P variants in PTPN11 (three), LZTR1 (three), and MRAS (one) were identified in six SPGF cases (including 4/31 GEMINI cases with CR) and one normozoospermic man. Undiagnosed RASopathies were detected in total for 17 ESTAND and GEMINI subjects, 15 SPGF patients (10 with CR), and two fertile men. Affected RASopathy genes showed high expression in spermatogenic and testicular somatic cells. In conclusion, congenital defects in the Ras/MAPK pathway genes represent a new congenital etiology of syndromic male infertility. Undiagnosed RASopathies were especially enriched among patients with a history of cryptorchidism. Given the relationship between RASopathies and other conditions, infertile men found to have this molecular diagnosis should be evaluated for known RASopathy-linked health concerns, including specific rare malignancies.

Genetic mutation of Cep76 results in male infertility due to abnormal sperm tail composition.

The transition zone is a specialised gate at the base of cilia/flagella, which separates the ciliary compartment from the cytoplasm and strictly regulates protein entry. We identified a potential new regulator of the male germ cell transition zone, CEP76. We demonstrated that CEP76 was involved in the selective entry and incorporation of key proteins required for sperm function and fertility into the ciliary compartment and ultimately the sperm tail. In the mutant, sperm tails were shorter and immotile as a consequence of deficits in essential sperm motility proteins including DNAH2 and AKAP4, which accumulated at the sperm neck in the mutant. Severe annulus, fibrous sheath, and outer dense fibre abnormalities were also detected in sperm lacking CEP76. Finally, we identified that CEP76 dictates annulus positioning and structure. This study suggests CEP76 as a male germ cell transition zone protein and adds further evidence to the hypothesis that the spermatid transition zone and annulus are part of the same functional structure.

The human infertility single-cell testis atlas (HISTA): an interactive molecular scRNA-Seq reference of the human testis.

BACKGROUND: Single-cell RNA-seq (scRNA-Seq) has been widely adopted to study gene expression of the human testis. Several datasets of scRNA-Seq from human testis have been generated from different groups processed with different informatics pipelines. An integrated atlas of scRNA-Seq expression constructed from multiple donors, developmental ages, and fertility states would be widely useful for the testis research community. OBJECTIVE: To describe the generation and use of the human infertility single-cell testis atlas (HISTA), an interactive web tool for understanding human spermatogenesis through scRNA-Seq analysis. METHODS: We obtained scRNA-Seq datasets derived from 12 donors, including healthy adult controls, juveniles, and several infertility cases, and reprocessed these data using methods to remove batch effects. Using Shiny, an open-source environment for data visualization, we created numerous interactive tools for exploring the data, some of which support simple statistical hypothesis testing. We used the resulting HISTA browser and its underlying data to demonstrate HISTA's value for testis researchers. RESULTS: A primary application of HISTA is to search by a single gene or a set of genes; thus, we present various analyses that quantify and visualize gene expression across the testis cells and pathology. HISTA also contains machine-learning-derived gene modules ("components") that capture the entire transcriptional landscape of the testis tissue. We show how the use of these components can simplify the highly complex data in HISTA and assist with the interpretation of genes with unknown functions. Finally, we demonstrate the diverse ways HISTA can be used for new data analysis, including hypothesis testing. DISCUSSION AND CONCLUSIONS: HISTA is a research environment that can help scientists organize and understand the high-dimensional transcriptional landscape of the human testis. HISTA has already contributed to published testis research and can be updated as needed with input from the research community or downloaded and modified for individual needs.

Toward clinical exomes in diagnostics and management of male infertility.

Infertility, affecting ∼10% of men, is predominantly caused by primary spermatogenic failure (SPGF). We screened likely pathogenic and pathogenic (LP/P) variants in 638 candidate genes for male infertility in 521 individuals presenting idiopathic SPGF and 323 normozoospermic men in the ESTAND cohort. Molecular diagnosis was reached for 64 men with SPGF (12%), with findings in 39 genes (6%). The yield did not differ significantly between the subgroups with azoospermia (20/185, 11%), oligozoospermia (18/181, 10%), and primary cryptorchidism with SPGF (26/155, 17%). Notably, 19 of 64 LP/P variants (30%) identified in 28 subjects represented recurrent findings in this study and/or with other male infertility cohorts. NR5A1 was the most frequently affected gene, with seven LP/P variants in six SPGF-affected men and two normozoospermic men. The link to SPGF was validated for recently proposed candidate genes ACTRT1, ASZ1, GLUD2, GREB1L, LEO1, RBM5, ROS1, and TGIF2LY. Heterozygous truncating variants in BNC1, reported in female infertility, emerged as plausible causes of severe oligozoospermia. Data suggested that several infertile men may present congenital conditions with less pronounced or pleiotropic phenotypes affecting the development and function of the reproductive system. Genes regulating the hypothalamic-pituitary-gonadal axis were affected in >30% of subjects with LP/P variants. Six individuals had more than one LP/P variant, including five with two findings from the gene panel. A 4-fold increased prevalence of cancer was observed in men with genetic infertility compared to the general male population (8% vs. 2%; p = 4.4 × 10-3). Expanding genetic testing in andrology will contribute to the multidisciplinary management of SPGF.

C19ORF84 connects piRNA and DNA methylation machineries to defend the mammalian germ line.

In the male mouse germ line, PIWI-interacting RNAs (piRNAs), bound by the PIWI protein MIWI2 (PIWIL4), guide DNA methylation of young active transposons through SPOCD1. However, the underlying mechanisms of SPOCD1-mediated piRNA-directed transposon methylation and whether this pathway functions to protect the human germ line remain unknown. We identified loss-of-function variants in human SPOCD1 that cause defective transposon silencing and male infertility. Through the analysis of these pathogenic alleles, we discovered that the uncharacterized protein C19ORF84 interacts with SPOCD1. DNMT3C, the DNA methyltransferase responsible for transposon methylation, associates with SPOCD1 and C19ORF84 in fetal gonocytes. Furthermore, C19ORF84 is essential for piRNA-directed DNA methylation and male mouse fertility. Finally, C19ORF84 mediates the in vivo association of SPOCD1 with the de novo methylation machinery. In summary, we have discovered a conserved role for the human piRNA pathway in transposon silencing and C19ORF84, an uncharacterized protein essential for orchestrating piRNA-directed DNA methylation.

Taking identity-by-descent analysis into the wild: Estimating realized relatedness in free-ranging macaques.

Biological relatedness is a key consideration in studies of behavior, population structure, and trait evolution. Except for parent-offspring dyads, pedigrees capture relatedness imperfectly. The number and length of DNA segments that are identical-by-descent (IBD) yield the most precise estimates of relatedness. Here, we leverage novel methods for estimating locus-specific IBD from low coverage whole genome resequencing data to demonstrate the feasibility and value of resolving fine-scaled gradients of relatedness in free-living animals. Using primarily 4-6× coverage data from a rhesus macaque (Macaca mulatta) population with available long-term pedigree data, we show that we can call the number and length of IBD segments across the genome with high accuracy even at 0.5× coverage. The resulting estimates demonstrate substantial variation in genetic relatedness within kin classes, leading to overlapping distributions between kin classes. They identify cryptic genetic relatives that are not represented in the pedigree and reveal elevated recombination rates in females relative to males, which allows us to discriminate maternal and paternal kin using genotype data alone. Our findings represent a breakthrough in the ability to understand the predictors and consequences of genetic relatedness in natural populations, contributing to our understanding of a fundamental component of population structure in the wild.

Improved phenotypic classification of male infertility to promote discovery of genetic causes.

An increasing number of genes are being described in the context of non-syndromic male infertility. Linking the underlying genetic causes of non-syndromic male infertility with clinical data from patients is important to establish new genotype-phenotype correlations. This process can be facilitated by using universal nomenclature, but no standardized vocabulary is available in the field of non-syndromic male infertility. The International Male Infertility Genomics Consortium aimed at filling this gap, providing a standardized vocabulary containing nomenclature based on the Human Phenotype Ontology (HPO). The "HPO tree" was substantially revised compared with the previous version and is based on the clinical work-up of infertile men, including physical examination and hormonal assessment. Some causes of male infertility can already be suspected based on the patient's clinical history, whereas in other instances, a testicular biopsy is needed for diagnosis. We assembled 49 HPO terms that are linked in a logical hierarchy and showed examples of morphological features of spermatozoa and testicular histology of infertile men with identified genetic diagnoses to describe the phenotypes. This work will help to record patients' phenotypes systematically and facilitate communication between geneticists and andrologists. Collaboration across institutions will improve the identification of patients with the same phenotypes, which will promote the discovery of novel genetic causes for non-syndromic male infertility.

Deleterious genetic changes in AGTPBP1 result in teratozoospermia with sperm head and flagella defects.

Approximately 10%-15% of couples worldwide are infertile, and male factors account for approximately half of these cases. Teratozoospermia is a major cause of male infertility. Although various mutations have been identified in teratozoospermia, these can vary among ethnic groups. In this study, we performed whole-exome sequencing to identify genetic changes potentially causative of teratozoospermia. Out of seven genes identified, one, ATP/GTP Binding Protein 1 (AGTPBP1), was characterized, and three missense changes were identified in two patients (Affected A: p.Glu423Asp and p.Pro631Leu; Affected B: p.Arg811His). In those two cases, severe sperm head and tail defects were observed. Moreover, AGTPBP1 localization showed a fragmented pattern compared to control participants, with specific localization in the neck and annulus regions. Using murine models, we found that AGTPBP1 is localized in the manchette structure, which is essential for sperm structure formation. Additionally, in Agtpbp1-null mice, we observed sperm head and tail defects similar to those in sperm from AGTPBP1-mutated cases, along with abnormal polyglutamylation tubulin and decreasing △-2 tubulin levels. In this study, we established a link between genetic changes in AGTPBP1 and human teratozoospermia for the first time and identified the role of AGTPBP1 in deglutamination, which is crucial for sperm formation.

Genome sequencing of Pakistani families with male infertility identifies deleterious genotypes in SPAG6, CCDC9, TKTL1, TUBA3C, and M1AP.

BACKGROUND: There are likely to be hundreds of monogenic forms of human male infertility. Whole genome sequencing (WGS) is the most efficient way to make progress in mapping the causative genetic variants, and ultimately improve clinical management of the disease in each patient. Recruitment of consanguineous families is an effective approach to ascertain the genetic forms of many diseases. OBJECTIVES: To apply WGS to large consanguineous families with likely hereditary male infertility and identify potential genetic cases. MATERIALS AND METHODS: We recruited seven large families with clinically diagnosed male infertility from rural Pakistan, including five with a history of consanguinity. We generated WGS data on 26 individuals (3-5 per family) and analyzed the resulting data with a computational pipeline to identify potentially causal single nucleotide variants, indels, and copy number variants. RESULTS: We identified plausible genetic causes in five of the seven families, including a homozygous 10 kb deletion of exon 2 in a well-established male infertility gene (M1AP), and biallelic missense substitutions (SPAG6, CCDC9, TUBA3C) and an in-frame hemizygous deletion (TKTL1) in genes with emerging relevance. DISCUSSION AND CONCLUSION: The rate of genetic findings using the current approach (71%) was much higher than what we recently achieved using whole-exome sequencing (WES) of unrelated singleton cases (20%). Furthermore, we identified a pathogenic single-exon deletion in M1AP that would be undetectable by WES. Screening more families with WGS, especially in underrepresented populations, will further reveal the types of variants underlying male infertility and accelerate the use of genetics in the patient management.

TAD evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function.

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find 14% of all human TAD boundaries to be shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons compared to species-specific boundaries. CRISPR-Cas9 knockouts of an ultraconserved boundary in a mouse model lead to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in the upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations and showcases the functional importance of TAD evolution.

AXDND1 is required to balance spermatogonial commitment and for sperm tail formation in mice and humans.

Dynein complexes are large, multi-unit assemblies involved in many biological processes including male fertility via their critical roles in protein transport and axoneme motility. Previously we identified a pathogenic variant in the dynein gene AXDND1 in an infertile man. Subsequently we identified an additional four potentially compound heterozygous variants of unknown significance in AXDND1 in two additional infertile men. We thus tested the role of AXDND1 in mammalian male fertility by generating a knockout mouse model. Axdnd1-/- males were sterile at all ages but could undergo one round of histologically complete spermatogenesis. Subsequently, a progressive imbalance of spermatogonial commitment to spermatogenesis over self-renewal occurred, ultimately leading to catastrophic germ cell loss, loss of blood-testis barrier patency and immune cell infiltration. Sperm produced during the first wave of spermatogenesis were immotile due to abnormal axoneme structure, including the presence of ectopic vesicles and abnormalities in outer dense fibres and microtubule doublet structures. Sperm output was additionally compromised by a severe spermiation defect and abnormal sperm individualisation. Collectively, our data highlight the essential roles of AXDND1 as a regulator of spermatogonial commitment to spermatogenesis and during the processes of spermiogenesis where it is essential for sperm tail development, release and motility.

A naturally occurring variant of MBD4 causes maternal germline hypermutation in primates.

As part of an ongoing genome sequencing project at the Oregon National Primate Research Center, we identified a rhesus macaque with a rare homozygous frameshift mutation in the gene methyl-CpG binding domain 4, DNA glycosylase (MBD4). MBD4 is responsible for the repair of C > T deamination mutations at CpG dinucleotides and has been linked to somatic hypermutation and cancer predisposition in humans. We show here that MBD4-associated hypermutation also affects the germline: The six offspring of the MBD4-null dam have a fourfold to sixfold increase in de novo mutation burden. This excess burden was predominantly C > T mutations at CpG dinucleotides consistent with MBD4 loss of function in the dam. There was also a significant excess of C > T at CpA sites, indicating an important, unappreciated role for MBD4 to repair deamination in CpA contexts. The MBD4-null dam developed sustained eosinophilia later in life, but we saw no other signs of neoplastic processes associated with MBD4 loss of function in humans nor any obvious disease in the hypermutated offspring. This work provides the first evidence for a genetic factor causing hypermutation in the maternal germline of a mammal and adds to the very small list of naturally occurring variants known to modulate germline mutation rates in mammals.

In vivo versus in silico assessment of potentially pathogenic missense variants in human reproductive genes.

Infertility is a heterogeneous condition, with genetic causes thought to underlie a substantial fraction of cases. Genome sequencing is becoming increasingly important for genetic diagnosis of diseases including idiopathic infertility; however, most rare or minor alleles identified in patients are variants of uncertain significance (VUS). Interpreting the functional impacts of VUS is challenging but profoundly important for clinical management and genetic counseling. To determine the consequences of these variants in key fertility genes, we functionally evaluated 11 missense variants in the genes ANKRD31, BRDT, DMC1, EXO1, FKBP6, MCM9, M1AP, MEI1, MSH4 and SEPT12 by generating genome-edited mouse models. Nine variants were classified as deleterious by most functional prediction algorithms, and two disrupted a protein-protein interaction (PPI) in the yeast two hybrid (Y2H) assay. Though these genes are essential for normal meiosis or spermiogenesis in mice, only one variant, observed in the MCM9 gene of a male infertility patient, compromised fertility or gametogenesis in the mouse models. To explore the disconnect between predictions and outcomes, we compared pathogenicity calls of missense variants made by ten widely used algorithms to 1) those annotated in ClinVar and 2) those evaluated in mice. All the algorithms performed poorly in terms of predicting the effects of human missense variants modeled in mice. These studies emphasize caution in the genetic diagnoses of infertile patients based primarily on pathogenicity prediction algorithms and emphasize the need for alternative and efficient in vitro or in vivo functional validation models for more effective and accurate VUS description to either pathogenic or benign categories.

Consensus label propagation with graph convolutional networks for single-cell RNA sequencing cell type annotation.

MOTIVATION: Single-cell RNA sequencing (scRNA-seq) data, annotated by cell type, is useful in a variety of downstream biological applications, such as profiling gene expression at the single-cell level. However, manually assigning these annotations with known marker genes is both time-consuming and subjective. RESULTS: We present a Graph Convolutional Network (GCN)-based approach to automate the annotation process. Our process builds upon existing labeling approaches, using state-of-the-art tools to find cells with highly confident label assignments through consensus and spreading these confident labels with a semi-supervised GCN. Using simulated data and two scRNA-seq datasets from different tissues, we show that our method improves accuracy over a simple consensus algorithm and the average of the underlying tools. We also compare our method to a nonparametric neighbor majority approach, showing comparable results. We then demonstrate that our GCN method allows for feature interpretation, identifying important genes for cell type classification. We present our completed pipeline, written in PyTorch, as an end-to-end tool for automating and interpreting the classification of scRNA-seq data. AVAILABILITY AND IMPLEMENTATION: Our code for conducting the experiments in this paper and using our model is available at https://github.com/lewinsohndp/scSHARP.

Rhesus macaque fetal and placental growth demographics: A resource for laboratory animal researchers.

Rhesus macaques (Macaca mulatta) are amongst the most common nonhuman primate species used in biomedical research. These animals provide a precious resource for translational studies and opportunities to maximize rhesus data use are encouraged. Here we compile data produced from 10 years of investigator-driven pregnancy studies conducted at the Oregon National Primate Research Center (ONPRC). All pregnancies were generated within the consistent and reproducible protocols of the ONPRC time-mated breeding program. The data included are from control animals who did not experience in utero perturbations or experimental manipulations. A total of 86 pregnant rhesus macaques were delivered by cesarean section over a range of gestational days (G) 50 to G159 (where term is G165 ± 10 days in the rhesus macaque), with subsequent immediate tissue harvesting following standardized protocols. Fetal and placental growth measures, and all major organ weights are reported. All data are presented relative to gestational age for the entire cohort and in addition, data are stratified by fetal sex. The outcome is a large reference resource for use by laboratory animal researchers in future comparative fetal development studies.

Genomic study of TEX15 variants: prevalence and allelic heterogeneity in men with spermatogenic failure.

Introduction: Human spermatogenesis is a highly intricate process that requires the input of thousands of testis-specific genes. Defects in any of them at any stage of the process can have detrimental effects on sperm production and/or viability. In particular, the function of many meiotic proteins encoded by germ cell specific genes is critical for maturation of haploid spermatids and viable spermatozoa, necessary for fertilization, and is also extremely sensitive to even the slightest change in coding DNA. Methods: Here, using whole exome and genome approaches, we identified and reported novel, clinically significant variants in testis-expressed gene 15 (TEX15), in unrelated men with spermatogenic failure (SPGF). Results: TEX15 mediates double strand break repair during meiosis. Recessive loss-of-function (LOF) TEX15 mutations are associated with SPGF in humans and knockout male mice are infertile. We expand earlier reports documenting heterogeneous allelic pathogenic TEX15 variants that cause a range of SPGF phenotypes from oligozoospermia (low sperm) to nonobstructive azoospermia (no sperm) with meiotic arrest and report the prevalence of 0.6% of TEX15 variants in our patient cohort. Among identified possible LOF variants, one homozygous missense substitution c.6835G>A (p.Ala2279Thr) co-segregated with cryptozoospermia in a family with SPGF. Additionally, we observed numerous cases of inferred in trans compound heterozygous variants in TEX15 among unrelated individuals with varying degrees of SPGF. Variants included splice site, insertions/deletions (indels), and missense substitutions, many of which resulted in LOF effects (i.e., frameshift, premature stop, alternative splicing, or potentially altered posttranslational modification sites). Conclusion: In conclusion, we performed an extensive genomic study of familial and sporadic SPGF and identified potentially damaging TEX15 variants in 7 of 1097 individuals of our combined cohorts. We hypothesize that SPGF phenotype severity is dictated by individual TEX15 variant's impact on structure and function. Resultant LOFs likely have deleterious effects on crossover/recombination in meiosis. Our findings support the notion of increased gene variant frequency in SPGF and its genetic and allelic heterogeneity as it relates to complex disease such as male infertility.

The origins and functional effects of postzygotic mutations throughout the human life span.

Postzygotic mutations (PZMs) begin to accrue in the human genome immediately after fertilization, but how and when PZMs affect development and lifetime health remain unclear. To study the origins and functional consequences of PZMs, we generated a multitissue atlas of PZMs spanning 54 tissue and cell types from 948 donors. Nearly half the variation in mutation burden among tissue samples can be explained by measured technical and biological effects, and 9% can be attributed to donor-specific effects. Through phylogenetic reconstruction of PZMs, we found that their type and predicted functional impact vary during prenatal development, across tissues, and through the germ cell life cycle. Thus, methods for interpreting effects across the body and the life span are needed to fully understand the consequences of genetic variants.

DDX3Y is likely the key spermatogenic factor in the AZFa region that contributes to human non-obstructive azoospermia.

Non-obstructive azoospermia, the absence of sperm in the ejaculate due to disturbed spermatogenesis, represents the most severe form of male infertility. De novo microdeletions of the Y-chromosomal AZFa region are one of few well-established genetic causes for NOA and are routinely analysed in the diagnostic workup of affected men. So far, it is unclear which of the three genes located in the AZFa chromosomal region is indispensible for germ cell maturation. Here we present four different likely pathogenic loss-of-function variants in the AZFa gene DDX3Y identified by analysing exome sequencing data of more than 1,600 infertile men. Three of the patients underwent testicular sperm extraction and revealed the typical AZFa testicular Sertoli cell-only phenotype. One of the variants was proven to be de novo. Consequently, DDX3Y represents the AZFa key spermatogenic factor and screening for variants in DDX3Y should be included in the diagnostic workflow.

TAD Evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function.

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species, and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find that only 14% of all human TAD boundaries are shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons, compared to species-specific boundaries. CRISPR-Cas9 knockouts of two ultraconserved boundaries in mouse models leads to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations, and showcase the functional importance of TAD evolution.

Diverse monogenic subforms of human spermatogenic failure.

Non-obstructive azoospermia (NOA) is the most severe form of male infertility and typically incurable. Defining the genetic basis of NOA has proven challenging, and the most advanced classification of NOA subforms is not based on genetics, but simple description of testis histology. In this study, we exome-sequenced over 1000 clinically diagnosed NOA cases and identified a plausible recessive Mendelian cause in 20%. We find further support for 21 genes in a 2-stage burden test with 2072 cases and 11,587 fertile controls. The disrupted genes are primarily on the autosomes, enriched for undescribed human "knockouts", and, for the most part, have yet to be linked to a Mendelian trait. Integration with single-cell RNA sequencing data shows that azoospermia genes can be grouped into molecular subforms with synchronized expression patterns, and analogs of these subforms exist in mice. This analysis framework identifies groups of genes with known roles in spermatogenesis but also reveals unrecognized subforms, such as a set of genes expressed across mitotic divisions of differentiating spermatogonia. Our findings highlight NOA as an understudied Mendelian disorder and provide a conceptual structure for organizing the complex genetics of male infertility, which may provide a rational basis for disease classification.